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Starting  from  published  charge/discharge  curves  and  X-ray  data  on  Pmmn-LiFe02  and  LiFe508  as  cath¬ 
ode  materials  vs.  Li  anode,  a  scheme  of  electrochemical  reactions  is  proposed  to  explain  the  unclear 
electrode  functionality  of  the  ‘corrugated  layer’  LiFe02  phase.  The  scheme  was  validated  by  quantum- 
mechanical  calculations  (CRYSTAL09  code,  hybrid  B3LYP  Hamiltonian)  on  a  number  of  structural 
models  for  Lii_xFe02,  LiFesOs,  and  Li3Fe508.  Magnetic  interactions  were  taken  into  account,  find¬ 
ing  antiferromagnetic  (Lii_xFe02)  and  ferrimagnetic  (LiFe508  and  Li3Fe508)  orderings  as  stable  states. 
At  variance  with  spinel-like  LiFe508,  Li3Fe508  displays  a  rocksalt-type  superstructure.  The  computed 
energies  for  reactions  (I)  4LiFe02  -+  4Li0.75FeO2  +  Li,  (II)  4Li0.75FeO2 +Li^  4/5LiFe508 +8/5Li20,  and  (III) 
1  / 2LiFe508  +  Li  4*  1  / 2Li3Fes08  are  4.44,  -3.62,  and  -2.1 0  eV,  respectively.  Such  values  compare  satisfacto¬ 
rily  with  the  average  charge/discharge  voltages  observed  for  positive  electrodes  made  up  of  Pmmn-LiFe02 
and  ofLiFe508. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  most  rechargeable  Li-ion  batteries  of  commercial  use,  the 
cathode  is  made  by  lithium  cobalt  oxide  against  an  anode  of 
graphite.  This  actually  works  satisfactorily,  but  two  inconveniences 
are  well  known:  cobalt  is  toxic  and  expensive,  so  as  to  raise 
problems  with  the  disposal  of  exhausted  devices  and  with  their 
production  cost.  Many  alternative  cathode  materials  have  been 
searched  for  in  the  last  two  decades  [1,2].  Those  based  on  lithium 
iron  oxides  are  particularly  attractive,  because  of  the  low  cost  and 
absence  of  toxicity  of  iron.  The  simplest  compound  of  this  family  is 
LiFe02,  which  indeed  shows  a  quite  complex  polymorphism. 

At  least  nine  phases  of  LiFe02  are  reported  in  the  literature 
[3].  The  first  group  is  formed  by  the  cubic  disordered  rocksalt- 
type  a  phase,  and  by  its  p  (tetragonal),  p'  (monoclinic),  and  7 
(tetragonal)  superstructures  [4,5].  Then  a  rhombohedral  ordered 
rocksalt  superstructure  analogous  to  LiCo02  (‘layer  structure’)  [6], 
a  hollandite-type  tetragonal  phase  [7],  and  three  orthorhombic 
phases  isostructural  with  LiMn02  [8],  with  p-NaFe02  [9],  and  with 
ramsdellite  [10]  were  described.  Yet  complete  structure  deter¬ 
minations  of  many  of  these  phases  are  missing.  Most  of  these 
polymorphs  were  tested  as  electrode  materials  against  a  lithium 
metal  anode.  The  claimed  performances  are  variable:  some  phases 
seem  not  to  work  (p'  and  7),  and  for  others  there  are  often  con¬ 
tradictory  reports.  Favourable  results  are  given  for  the  a  [11-13], 
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p-NaFe02-type  [9]  and  ramsdellite-type  [10,14]  modifications. 
The  best  reports  seem  to  be  those  concerning  the  so  called 
‘corrugated  layer’  polymorph,  with  orthorhombic  Pmmn  struc¬ 
ture  isomorphous  to  that  of  the  corresponding  phase  of  LiMn02 
(Fig.  1). 

Charge-discharge  curves  of  a  Pmmn-LiFe02  cathode  against 
Li  anode  were  determined  by  several  authors  [8,10,14-17].  In 
all  cases,  the  first  charge  cycle  shows  a  wide  voltage  plateau  at 
about  4.2  V,  suggesting  a  two-phase  character  of  the  electrode 
material  during  the  first  Li  extraction  cycle.  Then  the  following 
charge-discharge  cycles  take  a  more  regular  aspect,  with  the  volt¬ 
age  changing  in  the  1.5-3. 5  range  with  good  reversibility.  The 
reported  capacity  varies  from  about  1 00  mA  h  g-1  (0.3  0.4  e  f.u.-1 ) 

[8]  to  nearly  140mAhg-1  (0.5ef.u._1)  [17].  Different  interpreta¬ 
tions  are  given,  however,  on  the  basis  of  XRPD  data  collected  on  the 
electrode  material.  In  the  earlier  paper  (ex-situ  measurements),  an 
amorphization  process  of  the  cathode  during  the  first  charge  was 
suggested  [8].  On  the  other  hand,  in  situ  data  [17]  would  indicate 
the  formation  of  spinel-like  LiFesOg  in  the  first  discharge  cycle.  This 
interpretation  is  also  supported  by  magnetic  measurements,  show¬ 
ing  the  appearance  of  a  spontaneous  magnetic  moment  in  the  first 
discharge,  consistent  with  the  ferrimagnetic  character  of  LiFe508 
[17].  Further,  electrochemical  results  on  LiFesOs  demonstrated 
its  ability  to  intercalate/deintercalate  up  to  two  Li  atoms  f.u.-1 
reversibly  [18,19].  By  analysis  of  the  XRPD  intensities  collected 
after  discharge  [18],  the  formed  Li3Fe50g  phase  was  shown  to 
have  a  rocksalt-type  disordered  structure,  with  the  same  lattice 
as  spinel  but  with  a  different  cation  distribution.  This  feature 
could  not  appear  in  the  corresponding  electron  diffraction  pattern 
[19]. 
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Fig.  1.  Crystal  structure  of  Pmmn-LiFe02  [8].  The  Fe06  coordination  octahedra  and 
the  lithium  corrugated  layers  are  emphasized.  Li  and  O  atoms  are  dark  grey  (blue  in 
the  web  version  of  this  article)  and  white  balls,  respectively. 


On  the  basis  of  the  latter  experimental  data,  we  propose  the  fol¬ 
lowing  three-step  scheme  of  electrochemical  reactions  to  explain 
the  functionality  of  the  corrugated  layer  LiFe02  cathode. 

-  First  charge:  ~4.2V 

(1  /x)LiFe02  -+  (l/x)Lii_xFe02  +  Li  (I) 

An  x  <  1  \2  fraction  of  Li  is  extracted  from  lithium  iron  oxide,  with 
partial  oxidation  of  iron  and  formation  of  vacancies  in  the  Pmmn- 
LiFe02  structure. 

-  First  discharge:  4-2  V 

(l/x)Lii_xFe02  +  Li  -»  (l/5x)LiFe508  +  (2/5x)Li20  (II) 

Instead  of  reinserting  lithium,  the  Li-poor  iron  oxide  reacts 
with  lithium  metal  decomposing  into  the  spinel-like  phase  plus 
lithium  oxide,  with  reduction  of  Fe  to  the  original  +3  state.  Li20 
is  assumed  to  remain  in  the  cathode  as  electrochemically  inert 
phase,  without  taking  part  in  the  following  reactions.  Its  pres¬ 
ence  would  be  hard  to  detect  in  the  X-ray  pattern,  because  of 
peak  overlap  and  poor  scattering  power  of  Li. 

-  Subsequent  discharge/charge  cycles:  3.5-1. 5  V 

(l/2)LiFe508  +  Li  (l/2)Li3Fe508  (III) 

Lithium  is  intercalated/deintercalated  in  and  out  of  the  structure 
of  LiFe508,  with  the  oxidation  state  of  iron  switching  between  +3 
and  +2.6. 

Indeed,  a  final  short  plateau  at  about  4.5  V  is  also  observed  at 
the  end  of  the  charge  cycles  [17].  Although  this  feature  resembles 
the  long  continuous  plateau  in  the  first  charge,  which  was  ascribed 
to  reaction  (I),  delithiation  extends  here  well  beyond  the  previous 
capacity  range.  Further,  diffraction  and  magnetic  results  [  1 7  ]  do  not 
support  the  presence  of  LiFe02  in  the  electrode  after  cycling,  so  that 
a  repeat  of  reaction  (I)  is  not  likely.  Oxidation  of  LiFes08  to  some 
partially  delithiated  Li!_xFe508  species  would  be  probably  respon¬ 
sible  for  this  secondary  feature,  and  also  for  the  similar  charging  to 
4.5  Vreported  for  LiFes08  itself  [19].  Flowever,  this  process  will  not 
be  further  considered  in  the  present  work. 

In  order  to  validate  the  above  three-step  scheme,  and  thus  to 
explain  the  behaviour  of  Pmmn-LiFe02  as  cathode  material  on  a 
sound  basis,  a  first-principles  quantum-mechanical  study  of  reac¬ 
tions  (I)— (III)  and  of  all  phases  involved  was  undertaken.  Theoretical 
tools  proved  to  be  quite  effective  at  modelling  several  electrochem¬ 
ical  systems  [20-24],  thus  giving  an  important  contribution  to  their 
understanding  at  an  atomistic  level.  In  view  of  the  complexity  of  the 


present  case,  free  energy  changes  will  be  approximated  by  energy 
changes:  an  estimate  of  configurational  entropy  by  the  ideal  solu¬ 
tion  model  shows  that  its  contribution  is  comparatively  small  at 
room  temperature,  as  discussed  in  the  following  sections.  The  main 
emphasis  is  on  developing  suitable  structural  models  for  defective 
Lii_xFe02  and  for  Li3Fe508,  accounting  also  for  magnetic  ordering 
which  gives  a  substantial  energy  contribution  to  these  systems.  The 
final  aim  is  to  compute  the  reactions  energies  and  the  correspond¬ 
ing  average  electrochemical  potentials,  to  be  compared  with  the 
experimental  charge/discharge  voltages. 


2.  Computational  method 

Quantum-mechanical  calculations  of  the  ground-state  total 
crystal  energy  were  performed  by  the  computer  code  CRYSTAL09 
[25],  based  on  the  periodic  LCAO  (Linear  Combination  of  Atomic 
Orbitals)  approach.  A  hybrid  B3LYP  functional,  i.e.,  a  balanced 
mixture  of  the  DFT  ( Density-Functional-Theory  )-LYP  non-local  cor¬ 
relation  [26]  with  the  DFT  Becke’s  [27]  and  the  Flartree-Fock 
exchange,  was  employed.  All-electron  basis  sets  of  Gaussian-type 
functions  were  adopted  for  the  radial  parts  of  atomic  orbitals, 
according  to  the  schemes  8(s)641 1 1  (sp)41 1  (d)G,  8(s)41 1  (sp)l (d)G, 
and  5(s)ll(sp)G  for  Fe,  O  and  Li  atoms,  respectively,  in  all 
oxide  phases.  For  Li  metal,  the  richer  5(s)lll(sp)G  basis  set  was 
employed.  The  SCF  (Self-Consistent-Field)  equations  were  solved 
in  spin-polarized  form,  so  as  to  account  for  magnetic  ordering. 

The  reciprocal  space  was  sampled  according  to  a  regular  sublat¬ 
tice  of  volume  about  2x1 0-4  A-3  per  point.  For  instance,  in  the  case 
of  Pmmn-LiFe02  (a  =  4.061,  b  =  2.962  and  c  =  6.032  A  from  experi¬ 
ment  [  8  ] )  the  Monkhorst  grid  was  defined  by  4, 6,  and  3  points  along 
the  three  reciprocal  axes,  in  order  to  achieve  that  point  density. 
The  five  tolerances  related  to  cut-off  limits  for  series  summation 
were  set  to  10-7,  10-7,  10-7,  10-7,  and  10-14.  Convergence  was 
also  controlled  by  an  A E  threshold  of  10_9hartree  per  primitive 
unit  cell  in  the  SCF  cycles.  The  integration  of  the  DFT  functionals 
was  performed  by  use  of  the  LGRID  accuracy  conditions.  In  order 
to  accelerate  the  SCF  convergence,  the  technique  of  level  shifter 
was  used,  enhancing  the  energy  difference  between  highest  occu¬ 
pied  and  lowest  empty  states  in  the  first  cycles.  Atomic  coordinates 
and  unit-cell  constants  were  always  optimized  together  by  calcu¬ 
lation  of  analytical  gradients  and  subsequent  conjugate  gradients 
algorithm  (FULLOPTG  option). 

The  free  energy  changes  were  approximated  by  energy  changes 
at  0  K.  Thermal  and  configurational  entropy  effects  were  neglected. 
However,  an  estimate  of  configurational  entropy  according  to 
the  ideal  solution  model  shows  that,  for  reaction  (I),  at  room 
temperature  -TArS=  -0.058  and  -0.096 eV  for  x  =  l/4  and  1/16, 
respectively.  Such  values  are  less  than  2%  of  the  corresponding  A rE 
results  (Table  1). 


3.  Results  and  discussion 

3.1.  Vacancy  patterns  in  LU_xFe02  and  reaction  (I) 

In  order  to  do  calculations  for  the  first  reaction,  we  need  struc¬ 
tural  models  of  the  Li-deficient  Lii_xFe02  phase.  By  employing 
supercells  of  the  basic  stoichiometric  structure  (Fig.  1 ),  a  number  of 
models  were  built  up  for  each  of  the  compositions  x  =  1  /4, 1/8, 1/12, 
1/16.  For  each  of  them,  the  structure  was  optimized  by  full  relax¬ 
ation  obtaining  the  corresponding  total  energy  (Table  1 ).  The  local 
structural  distortions  in  the  vacancy  neighbourhood  can  be  clearly 
seen  in  Fig.  2,  for  three  of  the  four  vacancy  compositions.  How¬ 
ever,  the  changes  of  lattice  parameters  are  comparatively  small 
with  respect  to  pure  Pmmn-LiFe02.  The  simpler  ferromagnetic  (FM) 
ordering  was  assumed  in  all  cases,  and  the  antiferromagnetic  ( AFM) 
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Table  1 

Computed  energies  and  lattice  constants  of  the  least-energy  structural  models  of  Lii  _xFe02  for  different  vacancy  fractions  x.  Percent  deviations  from  the  experimental  cell 
edges  of  Pmmn-U¥e02  [8]  are  given  in  parentheses.  Ferromagnetic  (FM)  and  antiferromagnetic  (AFM)  ordering  is  indicated.  Z  is  the  number  of  formula  units  (f.u.)  per  unit-cell, 
and  A rE  refers  to  reaction  (I). 


X 

0 

1/16 

1/12 

1/8 

1/4 

FM  supercell 

do  x  bo  x  Co 

2  d0  x  4b0  x  c0 

2do  x  3bo  x  Co 

2do  x  2bo  x  Co 

2do  x  2bo  x  Co 

Space  group 

Pmmn 

Pmm2 

Pmm2 

Pmm2 

Cmm2 

Z 

2 

16 

12 

8 

8 

a  (A) 

4.101  (+1.0%) 

8.214 

8.219 

8.237 

8.291 

b 

2.975  (+0.4%) 

11.886 

8.895 

5.916 

5.882 

c 

6.124  (+1.5%) 

6.139 

6.143 

6.132 

6.143 

£+1419  (hartreef.ur1) 

-2.79436 

-2.31278 

-2.15224 

-1.83316 

-0.87257 

Ar£  (eV) 

4.825 

4.830 

4.399 

4.332 

A(x)  (eV) 

0 

0.005 

-0.426 

-0.493 

AFM  supercell 

2ao  x  h0  x  Co 

2do  x  2bo  x  Co 

S.G. 

P21/m 

Pm 

Z 

4 

8 

a  (A) 

8.104 

8.163 

b 

2.982 

5.876 

c 

6.163 

6.233 

P (deg) 

90.21 

89.97 

£+1419  (hartreef.u.-1) 

-2.79730 

-0.87453 

Ar£  (eV) 

4.439 

£  (AFM)  -  £  (FM)  (eV  f.u.-1 ) 

-0.080 

-0.053 

Fig.  2.  Lithium  vacancy  (□)  arrangements  in  the  (001)  corrugated  layers  of  the  least-energy  optimized  structures,  for  three  compositions  (x  =  1/4, 1/8,  and  1/12)  of  Lii  _xFeC>2. 
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one  was  considered  only  for  the  two  compositional  end-members. 
We  shall  discuss  first  the  FM  results,  where  the  relevant  physics  of 
the  system  is  already  accounted  for. 

The  most  significant  quantity  reported  in  Table  1  is: 

A  rE(x)  =  (l/x)[£(Lii_xFe02)-£(LiFe02)]  +  £(Li)  =  (l/x)AE  +  E(Li) 

This  is  the  reaction  (I)  energy,  i.e.  the  energy  necessary  to  extract 
Li  atoms  with  x  density  from  the  Pmmn-LiFe02  structure,  and  to 
insert  them  into  BCC  lithium  metal,  per  Li  atom.  It  also  corre¬ 
sponds  to  the  formation  energy  of  a  Li  vacancy  with  density  x 
in  Pmmn-LiFe02,  assuming  as  energy  zero  that  of  lithium  metal. 
The  calculated  value  of  E(Li)  was  -7.527956  hartree  atom-1,  with 
a  =  3.420  A  as  optimized  cell  edge.  By  normalizing  reaction  (I)  to  one 
Li  atom  one  electron  is  exchanged  in  the  Redox  process,  so  that  the 
A rE  value  measured  in  eV  is  numerically  equal  to  the  opposite  of 
the  electrochemical  potential  E  measured  in  Volts. 

The  dependence  of  A rE  on  x  indicates  the  effect  of  vacancy- 
vacancy  interactions  on  their  formation  energy.  The  structural  con¬ 
figurations  of  the  two  diluted  vacancy  cases  (x  =  l/16  and  1/12) 
approximate  the  isolated  vacancy  situation  (cf.  Fig.  2),  where  inter¬ 
actions  are  absent.  The  corresponding  4.83  eV  value  can  thus  be 
considered  as  an  estimate  of  the  ideal  single  vacancy  formation 
energy.  Deviations  from  this  value,  expressed  as  A(x)  in  Table  1 ,  are 
a  measure  of  vacancy-vacancy  interactions.  Both  vacancy-richer 
compositions  (x  =  l/8  and  1/4)  show  negative  values  of  -0.4  to 
-0.5  eV,  indicating  a  stabilization  interaction.  This  is  due  to  the 
peculiar  arrangement  of  [0 1  0]  rows  of  alternating  Li  atoms  and 
vacancies,  which  is  observed  in  both  cases  (Fig.  2).  It  thus  appears 
that  a  pair  of  vacancies  in  a  [0 1  0]  row  is  more  stable  when  sepa¬ 
rated  by  a  single  Li  atom,  than  by  two  or  more  of  them.  This  was 
confirmed  by  building  up  structural  models  of  the  vacancy-rich 
compositions  where  such  an  arrangement  was  absent,  and  their 
energies  proved  to  be  similar  to  those  of  the  vacancy-poor  cases. 
For  instance,  for  x  =  1/8  a  Cmm2  model  with  2ao  x  4 b0  x  Co  super¬ 
cell  has  vacancies  separated  by  three  Li  atoms  along  the  row:  the 
computed  energy  was  4.95  eV,  with  a  slight  positive  deviation  with 
respect  to  isolated  vacancies.  On  the  other  hand,  the  configuration 
with  two  adjacent  vacancies  along  the  row  gave  a  large  destabiliza¬ 
tion  (x=  1/4,  Pmm2  a0  x  4 b0  x  c0  supercell,  ArE  =  5.23  eV).  We  can 
then  conclude  that,  for  phases  with  large  vacancy  density  (x  >  0.1 ), 
the  [0 1  0]  Li-D-Li-D-Li  arrangement  should  be  present  and  the 
vacancy  formation  energy  A rE  should  be  similar  to  the  4.33  eV 
value  obtained  for  x  =  0.25. 

Because  of  the  two-phase  character  of  the  first  charge  curve 
obtained  in  the  experimental  electrochemical  cell,  a  composition 
discontinuity  can  be  assumed  between  stoichiometric  Pmmn- 
LiFe02  and  the  Lii_xFe02  phase  with  the  largest  Li  deficiency 
(x  >  0.25).  Therefore,  the  corresponding  reaction  energy  of  4.33  eV 
is  the  appropriate  constant  value  for  the  whole  x  range.  This  energy 
compares  quite  well  with  the  experimental  4.2  V  voltage  of  the 
plateau  observed  in  the  first  charge  cycle. 

Pmmn-LiFe02  is  reported  to  be  antiferromagnetic  from  0  K  to 
room  temperature,  on  the  basis  of  incomplete  magnetic  data  [8,17]. 
Then  the  simplest  AFM  ordering  scheme  was  adopted  for  the  LiFe02 
structure:  the  a0  unit-cell  edge  was  doubled,  so  as  to  alternate  Fe3+t 
and  Fe3+|  ions  along  the  [10  0]  rows  of  corner-sharing  Fe06  octahe- 
dra,  where  Fet-O-Fej,  superexchange  interactions  can  take  place. 
The  same  scheme  holds  for  Li0j5FeO2,  too.  As  expected,  stabiliza¬ 
tion  with  respect  to  the  FM  arrangement  is  observed  for  both  cases 
(Table  1 );  however,  the  effect  is  smaller  for  the  Li-deficient  compo¬ 
sition.  This  means  that  superexchange  is  perturbed  somehow  by  Li 
vacancies.  As  a  result,  the  reaction  energy  increases  slightly  from 
4.33  to  4.44  eV;  this  value  should  be  the  preferred  one,  as  it  is  based 
on  the  most  stable  magnetic  configurations 


Table  2 

Splitting  of  Wyckoff  sites  of  the  spinel  structure  with  symmetry  lowering  down  to 
orthorhombic  C222i. 


Space  group 

Fd3m  ->  P4332 

Octahedral  sites 

16d 

12di  +4b 

16c 

12d2  +  4a 

Tetrahedral  sites 

8  a 

8ca 

8  b 

8cb 

Space  group 

P4332  ->  C222i 

Octahedral  sites 

12di 

8ci  +8c2  +4ai  +4bi 

4b 

4  a2  +  4  b2 

12  d2 

8c3  +  8C4  +  4a3  +  4  b3 

4a 

4a4  +  4b  4 

Tetrahedral  sites 

8ca 

8C5  +  8C6 

3.2.  Magnetic  structures  ofLiFesOg  and  Li3Fe508 

The  crystal  structure  of  LiFesOg  (space  group  P4332)  is  an 
ordered  version  of  inverse  spinel  (Fd3m),  with  Li  and  3/5  of  Fe  atoms 
in  octahedral  sites,  and  2/5  of  Fe  atoms  in  tetrahedral  sites  (Table  2) 
[28].  Magnetic  measurements  indicate  that  this  compound  is  fer- 
rimagnetic  (FFM)  at  room  temperature  [29-31]:  the  net  magnetic 
moment  arises  then  from  the  difference  between  octahedral  Fe3+ 
ions  with  majority  spin  and  tetrahedral  Fe3+  ions  with  minority  spin 
(Fig.  3).  Thus,  a  study  of  the  magnetic  structure  of  LiFesOg  should  be 
included  in  the  modelling  to  insure  appropriate  energy  results.  Both 
the  FM  and  FFM  magnetic  orderings  are  compatible  with  the  P4332 
space  group,  but  any  antiferromagnetic  arrangement  of  the  spins 
requires  that  the  symmetry  is  lowered  to  a  subgroup  of  P4332.  The 
simplest  AFM  model  can  be  built  up  within  the  C222!  symmetry 
with  a  ^2a0  x  ^2a0  x  a0  conventional  cell  (Table  2);  with  respect 
to  the  FFM  case,  1  /6  of  the  octahedral  Fe3+  ions  move  from  majority 
to  minority  spin,  so  as  to  make  the  net  magnetic  moment  vanish. 

In  Table  3,  the  ideal  fractional  coordinates  of  the  C222i  struc¬ 
ture  are  reported,  with  the  site  occupation  pattern  of  AFM-LiFesOg. 
The  site  patterns  of  all  three  magnetic  structures  (FM,  FFM  and 
AFM)  considered  for  LiFesOg,  in  the  corresponding  space  groups, 
are  given  in  Table  4.  Therein  also  the  results  of  the  least-energy 
structure  optimizations  appear  for  each  phase.  Very  strong  A E 
stabilization  energies  are  observed  for  both  the  FFM  and  AFM  order¬ 
ings  with  respect  to  the  FM  configuration;  the  values  are  about 
three  times  that  computed  for  AFM-LiFe02  (cf.  Table  1).  However, 


Fig.  3.  Picture  of  the  ferrimagnetic  (cubic  P4332)  structure  of  LiFe508.  Li  and  O  atoms 
are  denoted  as  dark  grey  (blue  in  the  web  version)  and  white  balls,  respectively. 
Fe  atoms  are  vertically-  (majority  spin)  and  horizontally-  (minority  spin)  hatched 
spheres. 
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Table  3 

Ideal  fractional  coordinates  of  the  sites  of  space  group  C222i  occupied  by  Fe  and  Li 
atoms  in  LiFe5Og  and  Li3Fe508.  All  sites  are  octahedral  but  those  denoted  by  aster¬ 
isks  (tetrahedral).  Antiferromagnetic  (AFM)  and  ferromagnetic  (FFM)  orderings  are 
specified. 


Site 

X 

y 

z 

AFM-LiFe5Og 

FFM-Li3Fe50g 

8d 

1/4 

1/8 

1/4 

Fet 

Fet 

8C2 

3/8 

1/4 

0 

Fet 

Fet 

8c3 

1/8 

1/4 

0 

Li 

8C4 

1/4 

3/8 

1/4 

Fet 

8c5* 

7/8 

7/8 

3/8 

Fet 

8c6* 

5/8 

1/8 

7/8 

Fet 

4a  i 

7/8 

0 

0 

Fet 

Fet 

4b  i 

1/2 

1/8 

1/4 

Fet 

Fet 

402 

3/8 

0 

0 

Li 

Li 

4  b2 

0 

1/8 

1/4 

Li 

Li 

4  a3 

5/8 

0 

0 

Fet 

4  b3 

0 

3/8 

1/4 

Li 

404 

1/8 

0 

0 

Fet 

4  b4 

1/2 

3/8 

1/4 

Li 

the  FFM  result  is  significantly  larger  in  modulus  (-0.262  against 
-0.21 2  eV  Fe  atom-1 )  than  the  AFM  one,  confirming  that  the  stable 
structure  of  LiFesOg  is  ferrimagnetic  as  shown  by  experiment.  The 
calculated  volume  of  the  FFM  phase  also  shows  the  best  agreement 
with  the  experimental  value  [18]  (+3.7%  deviation),  whereas  for  the 
other  phases  larger  positive  deviations  are  observed. 

In  order  to  accommodate  additional  Li  atoms  into  the  spinel-like 
LiFe508  structure,  the  4 a  and  12d2  empty  octahedral  sites  cannot 
be  utilized,  because  they  are  too  close  to  Eef  atoms  in  the  tetrahe¬ 
dral  8ca  position.  Similarly,  the  empty  tetrahedral  8 cb  site  is  useless 
because  near  to  Fet  in  \2d\  e  to  Li  in  4b.  Therefore,  Fet  atoms 
must  be  moved  from  the  tetrahedral  8ca  site  into  an  octahedral 
one,  so  as  to  insert  additional  lithium  into  the  other  empty  octa¬ 
hedral  positions.  The  structure  obtained  has  thus  all  metal  atoms 
in  octahedral  sites,  with  the  tetrahedral  ones  empty:  it  is  no  more 
a  spinel,  but  a  rocksalt-type  superstructure.  For  the  Li3Fe508  com¬ 
position,  the  only  possible  arrangement  within  the  original  P4332 
space  group  is  Fe  in  12di,  4a  and  4b,  and  Li  in  12d2.  Unfortunately, 
this  arrangement  is  not  consistent  with  the  Rietveld-refined  occu¬ 
pancies  of  0.5(Fe),  0.5(Li)  for  12d2  and  4a,  with  full  occupancies  of 


A 

z 


Fig.  4.  Ordered  (orthorhombic  C222i)  ferrimagnetic  structural  model  of  Li3Fe508. 
Black  balls  indicate  additional  intercalated  Li  atoms  with  respect  to  LiFe508. 

Fe  in  12di  and  Li  in  4 b  [18].  Further,  an  attempt  to  calculate  this 
P4332  structure  yielded  invariably  a  conductive  state  with  difficult 
SCF  convergence. 

The  simplest  ordered  arrangement  of  Li3Fe508  consistent  with 
the  average  disordered  P4332  structure  from  diffraction  results  can 
be  achieved  within  the  C222i  subgroup.  An  energy  screening  of  the 
nearly-degenerate  Fe/Li  distributions  among  sites  indicated  the  one 
reported  in  Tables  3  and  5  as  the  most  favoured  one  (Fig.  4).  Inser¬ 
tion  of  additional  Li  atoms  and  shift  of  iron  atoms  from  tetrahedral 
to  octahedral  sites,  on  passing  from  the  LiFe508  to  the  Li3Fe508 
compound,  are  clearly  shown  in  Fig.  5.  The  arrangement  of  Fe  atoms 
in  the  C222!  structure  of  Li3Fe508  is  compatible  with  both  FM  and 
FFM  orderings;  the  corresponding  schemes  are  reported  in  Table  5. 
The  results  of  the  least-energy  optimizations  show  that,  also  for 
Li3Fe508,  the  FFM  arrangement  is  the  most  stable  one.  Flowever, 
the  AE  =  E(FFM)  —  £(FM)  stabilization  energy  is  only  25%  of  that 
obtained  for  LiFe508  (cf.  Tables  4  and  5).  This  means  that  the  ferri¬ 
magnetic  spin  ordering  process  5Fet  -+  3Fet  +  2Fet  is  energetically 
more  favoured  when  the  chemical  environment  of  spin-up  (octa¬ 
hedral)  and  spin-down  (tetrahedral  sites)  is  different,  than  when  it 


Table  4 

Wyckoff  site  occupation  patterns,  calculated  lattice  constants  and  energies  of  the  magnetic  structures  of  LiFe5Og. 


FM  P4332 

FFM  P4332 

AFM  C222i 

Exp.  P4332  [18] 

Fet 

12di  +8ca 

12di 

8ci  +  8C2  +  4ai 

Fet 

8ca 

8c5  +  8c6  +  4bi 

Li 

4b 

4b 

4a2  +  4  b2 

a  (A) 

8.478 

8.415 

11.951 

8.331 

b 

11.911 

c 

8.414 

V  (A3  f.u.-1 ) 

152.34 

148.97 

149.71 

143.67 

E  +  6928  (hartree  f.u.-1 ) 

-0.01132 

-0.05948 

-0.05029 

AE  (eV  Fe  atom-1) 

0 

-0.262 

-0.212 

Table  5 

Wyckoff  site  occupation  patterns,  calculated  lattice  constants  and  energies  of  the  magnetic  structures  of  Li3Fe50g. 

FM  C222i 

FFM  C222x 

Exp.  P4332  [18] 

Fet 

8ci  +8C2  +4ai  +4bi 

8ci  +8c2  +  4ai  +4bi 

8c4  +  4a3  +  4a4 

Fet 

8c4  +  4a3  +  4a4 

Li 

4a2  +  4  b2 

4a2  +  4  b2 

8  c3  +  4  fo3  +  4  b4 

8c3  +  4  fo3  +  4  b4 

a  (A) 

12.023 

12.033 

8.389 

b 

12.011 

11.983 

c 

8.541 

8.526 

V  (A3  f.u.-1 ) 

154.17 

153.67 

147.60 

E  +  6943  (hartree  f.u.-1 ) 

-0.25704 

-0.26946 

AE  (eVFe  atom-1) 

0 

-0.068 
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Fig.  5.  (001)  atomic  layers  at  z  =  0  and  z=  1/4  of  (a)  LiFe508  and  (b)  Li3Fe508  in  the  C222\  reference  frame;  in  (a)  the  smaller  cubic  spinel  unit-cell  is  outlined,  too.  Cf. 
Figs.  2  and  4  for  notation. 


is  equal  (octahedral  sites  for  both  spins).  Consequences  ensue  for 
the  electrochemical  reaction  (III)  energy,  as  discussed  in  the  next 
section. 

3.3.  Energetics  of  reactions  (II)  and  (III) 

The  reaction  energy  of  transformation  (II)  was  computed  with 
reference  to  the  x  =  0.25  case,  according  to  the  expression  given  in 
Table  6,  on  the  basis  of  the  energies  previously  discussed  for  the 
phases  involved.  As  for  Li20  with  antifluorite  Fm3m  structure,  the 
E  =  -90.444576  hartree  f.u.-1  value  (with  a  =  4.545  A)  was  obtained 
and  employed.  The  values  of  ArE(II)  were  calculated  by  combin¬ 
ing  the  energies  of  the  different  magnetic  states  of  LiFe508  and 
Li0.75FeO2  (Tables  1  and  4),  to  show  the  effect  of  magnetic  order¬ 
ing  of  these  compounds  on  the  overall  reaction  energy,  and  they 
are  reported  in  Table  6.  The  most  reliable  result  (-3.615  eV)  cor¬ 


responds  to  the  FFM-LiFe5O8/AFM-Li0.75FeO2  case  involving  the 
most  stable  magnetic  structures,  in  agreement  also  with  experi¬ 
mental  observation.  The  computed  reaction  (II)  energy  fits  in  the 
4-2  V  range  of  the  electrochemical  potential  measured  in  the  first 
discharge  cycle,  though  it  is  slightly  larger  than  the  average  value. 

Similar  calculations  were  performed  for  reaction  (III),  on  the 
basis  of  the  total  energies  obtained  for  different  magnetic  arrange¬ 
ments  of  LiFe508  and  Li3Fe508  (Tables  4  and  5).  The  results  of 
A rE  are  reported  in  Table  6,  and  again  they  appear  to  be  sub¬ 
stantially  affected  by  the  kind  of  magnetic  ordering  assumed.  The 
FFM  vs.  FM  stabilization  is  much  smaller  for  Li3Fes08  than  for 
LiFesOs,  so  that  the  reaction  energy  is  significantly  reduced  (in 
absolute  value)  when  both  phases  are  FFM-ordered.  Since  this  is 
the  most  stable  arrangement,  the  corresponding  ArE=  -2.096  eV 
value  should  be  taken  as  correct  predicted  reaction  (III)  energy.  This 
has  to  be  compared  with  the  average  potential  of  the  experimen- 


Table  6 

Calculated  reaction  energy  (eVLi  atom-1)  for  the  two  electrochemical  reactions  supposed  to  occur  in  the  first  discharge  (II)  and  subsequent  charge-discharge  cycles  (III)  of 
the  Pmmn-LiFe02  cathode  vs.  a  Li  metal  anode.  Different  magnetic  orderings  of  LiFe508,  Li0.75FeO2  and  Li3Fe508  are  considered. 

(II)4Li0.75FeO2  +  Li^  (4/5)LiFe508  +  (8/5)Li20  ;  ArE  =  (4/5)£(LiFe508)  - 4£(Lio.75Fe02)  +  (8/5)£(Li20)  - £(Li) 

A rE  FM-FM  A rE  FFM-FM 

-2.780  -3.828 


(III)(l/2)LiFe508  +  Li*>(l/2)Li3Fe508  ;  ArE  =  (l/2)[E(LiFe508)-E(LiFe508)]  -E(Li) 

A  rE  FM-FM  A  rE  FFM -FFM 

-2.583  -2.096 


A  rE  AFM-AFM 
-3.414 


A  rE  FFM-AFM 
-3.615 
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tal  charge/discharge  curves  of  Pmmn- LiFe02  in  the  1.5-3  V  range 
[8,17],  and  to  those  of  LiFe508  in  the  1.2-3  V  [18]  or  1.7-3.5  V  [19]. 
The  corresponding  integral  averages  are  about  2.0, 1.6  and  2.2  V,  so 
that  our  predicted  energy  agrees  very  satisfactorily  with  the  more 
recent  experimental  results  [8,17,19]. 

4.  Conclusions 

A  multi-step  scheme  of  electrochemical  reactions  was  pro¬ 
posed  to  explain  charge/discharge  voltage  and  X-ray  diffraction 
data  of  Pmmn-LiFe02  as  positive  electrode  against  a  lithium  anode. 
In  particular,  formation  of  spinel-like  LiFesOg,  and  its  subse¬ 
quent  reversible  lithiation  to  a  rocksalt-type  phase,  are  the  key 
stages  of  the  process.  The  third  step  has  also  an  interest  of  its 
own,  as  it  concerns  the  operation  of  LiFe508  as  starting  cath¬ 
ode  material.  By  accurate  first-principles  calculations,  least-energy 
structural  models  were  determined  for  the  vacancy  distribution  in 
defective  Lii_*Fe02  and  for  Li3Fe508,  including  various  magnetic 
arrangements.  Antiferromagnetic  ordering  for  Lii_xFe02,  and  fer- 
rimagnetic  ordering  for  both  LiFe508  and  Li3Fe508  are  found  to 
characterize  the  corresponding  ground  states.  The  reaction  ener¬ 
gies  at  the  athermal  limit  were  then  computed  for  all  three  steps 
of  the  process,  obtaining  values  which  compare  very  favourably 
with  the  average  charge/discharge  voltages  reported  from  exper¬ 
iment.  This  supports  strongly  the  proposed  scheme  of  atomic 
mechanisms  for  the  ‘corrugated  layer’  lithium  iron  oxide  electrode. 
First-principles  calculations  are  confirmed  to  be  a  valuable  aid  for 
interpreting  and  understanding  electrode  functionality  in  batter¬ 
ies. 
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